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Abstract. Understanding the sources responsible for reionizing the Universe is a key goal of 
observational cosmology. A discrepancy has existed between the metagalactic hydrogen ioniza- 
tion rate, Thi, predicted by early hydrodynamical simulations of the Lyman-a forest if scaled 
to appropriate assumptions for the IGM temperature, when compared to values predicted from 
the proximity effect. We present new estimates for Thi in the redshift range 2 < z < 4 based 
on hydrodynamical simulations of the Lyman-a forest opacity. Within the current concordance 
cosmology, and assuming updated QSO emissivity rates, a substantial contribution to the UV 
background from young star-forming galaxies appears to be required over the entire redshift 
range. Our results are consistent with lower-end estimates from the proximity effect. It is also 
found that the errors on the ionization rate are dominated by uncertainties in the thermal state 
of the intergalactic medium and the r.m.s fluctuation amplitude at the Jeans scale. 
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1. Introduction 

The series of low column density (10 125 < Nm < 10 14 5 cm -2 ) Lya absorption features 
seen blue- ward of the Lya emission line in high redshift quasi-stellar object (QSO) spectra 
trace the distribution of neutral hydrogen in the intergalactic medium (IGM) . This forest 
of absorption lines provides a unique, unbiased probe of the thermal and ionization 
history of the IGM along the QSO line-of-sight. This has motivated many authors to use 
hydrodynamical simulations of structure formation, calibrated to reproduce the observed 
properties of the forest obtained from high resolution QSO spectra, to infer the amplitude 
of the metagalactic UV background (see the seminal paper by Rauch et al. 1997, and 
most recently Tytler et al. 2004). Consequently there are a wide range of estimates for 
the hydrogen ionization rate per atom, Thi , inferred from the Lya forest opacity using 
numerical simulations with different assumptions. This complicates the comparison with 
determinations of the ionization rate from estimates of the integrated emission from 
observed QSOs and/or galaxies (e.g. Haardt & Madau 1996, hereafter HM96; Madau, 
Haardt & Rees 1999) and estimates using the proximity effect (e.g. Scott et al. 2000) 
which both also have rather large uncertainties. Furthermore, Steidel et al. (2001), in a 
study of Lyman break galaxies at z ~ 3.4, suggested that the intensity of the ionizing 
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background may be a factor of a few larger than in the model of HM96 due to a large 
contribution from star-forming galaxies. 

We discuss new estimates of the metagalactic hydrogen ionization rate presented in 
Bolton et al. (2005), inferred from state-of-the-art hydrodynamical simulations of the 
Lya forest. These estimates are consistent with a substantial contribution from galax- 
ies to the amplitude of the UV background at the hydrogen ionization edge, and are 
in agreement with other observational estimates of the dimcnsionlcss quantity T_i2 — 

Wio- 1 ^- 1 . 

2. Numerical code and convergence 

We use a suite of 21 high resolution hydrodynamical simulations run using a new 
version of GADGET (Springel, Yoshida & White, 2001) to investigate in detail the de- 
pendence of Thi on physical and numerical parameters (for a detailed description of the 
simulations see Bolton et al. 2005). The cosmological parameters for our adopted fiducial 
model are: 

(fi m , n h h 2 , h, ct 8 ) = (0.26 ± 0.04, 0.024 ± 0.001, 0.72 ± 0.04, 0.85 ± 0.05), 

consistent with the results of Spergel et al. (2003), assuming a flat, vacuum energy 
dominated cosmological model. We assume n = 0.95 and a helium mass fraction of 
Y = 0.24. In addition, we adopt fiducial values for the thermal state of the gas at mean 
cosmic density; T = [11200, 17800, 12500] ± 5000 K at z = [2, 3, 4] and 7 = 1.3 ± 0.3 for 
the index of the gas effective equation of state (Hui & Gnedin 1997), based on the results 
of Schaye et al. (2000). We linearly rescale the opacity of the artificial Lya forest spectra 
we construct to match the central values for the Lya forest effective optical depth taken 
from the fitting formula of Schaye et al. (2003): r cff = [0.130±0.021, 0.362±0.036, 0.805± 
0.070]. We treat the amplitude of the UV background as a free parameter, so rescaling 
the neutral hydrogen opacity corresponds to linearly rescaling the hydrogen ionization 
rate by the same amount. 

We initially run 7 simulations of our fiducial model with differing box size and reso- 
lution to assess the degree of numerical convergence. We find simulations with less than 
400 3 gas particles within a 30/i _1 comoving Mpc box do not achieve adequate numerical 
convergence. We take our fiducial model to have 200 3 gas particles within a 15/i _1 co- 
moving Mpc box; we estimate the value of T_i2 we infer will be systematically high by 
around 8 per cent due to the combined error from box size and resolution. 

3. Scaling relations for T 12 from simulated absorption spectra 

When inferring T_i2 from hydrodynamical simulations and subsequently rescaling with 
different parameters, it is generally assumed that (e.g Rauch et al. 1997): 

r_i 2 OC O^3 T -0.7 O -0.5 (31) 

It is implicit in this relation that the density and velocity distribution, along with the 
effective equation of state of the Lya absorbers remain unchanged for different values 
of fib, h, fi m , T and r_i2. We test the validity of equation 13. II using 14 simulations in 
addition to our fiducial model, by varying the parameters T, fi m and as in each run. 
We also rescale to different values for the temperature-density relation index 7 and the 
effective optical depth in post-processing. We do not discuss the scaling with fib and h, 
since our analysis revealed that equation 13. II holds extremely well for these parameters. 
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Figure 1. Left: The dependence of T_i2 on the gas temperature at mean density To at the 
three different redshifts indicated on the plot. Solid curves show a least squares fit and the 
dotted curves shows the T -0,7 scaling due to the temperature dependence of the recombination 
coefficient. The filled symbols show the fiducial temperatures and their uncertainties. Right: The 
dependence of F_i2 on the index of the temperature density relation 7. The filled circle shows 
our fiducial value of 7 and its uncertainty. 



Table 1. The redshift dependent indices from our scaling relation of T_i2 with several 
cosmological and astrophysical parameters. 
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-0.90 
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-0.99 


-1.61 
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-0.59 
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-1.26 
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Figures and |21 and Table ^ summarise the main results of this study. We find that 
r_i2 scales around our fiducial model as: 

r_ 12 oc nlh 3 T x ^ z h X2(z) ^ (z) crt i{z) T2 {z \ (3.2) 

where we tabulate the indices of equation 13.21 in Table We stress that this scaling 
relation is likely to be somewhat model dependent and should not be applied to models 
with parameters very different from our fiducial model without further checks. However, 
it does provide a clear picture of the degeneracies which exist between T_i2 and other 
parameters within our simulations, and is independent of assumptions about the gas dis- 
tribution, effective equation of state and ionized gas fraction of the low density 1GM. In 
particular, the Lycv optical depth at fixed r.m.s. fluctuation amplitude as is more strongly 
dependent on f2 m than eauation l3. II suggests due to changes in the gas density distribu- 
tion. Thermal broadening also produces a deviation from the scaling of eauation l3.1l for 
temperature, although this change is much less dramatic. Consequently, we urge caution 
when using eauation l3.ll especially if comparing models with differing f2 m . We also find 
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Figure 2. Left: The dependence of the estimated r_i2 on fi m for the three different redshifts 
indicated on the plot. The solid curves are a least-square fit. The dotted lines show the f2m°' 5 
scaling. Filled points are obtained for a model with f2 m = 1.0 and the the same r.m.s. fluctuation 
amplitude at a scale of 30 kms" 1 as our fiducial model. The filled circle shows our fiducial value 
of f2 m = 0.26 and its uncertainty. Right: The dependence of F_i2 on og. The dotted curve shows 
the result of Tytler et al. at z = 1.9, assuming (F) b s = 0.882 ±0.01. The filled circle shows our 
fiducial value of erg and its uncertainty 

a strong dependence of T_i2 on the effective optical depth of the Lya forest which the 
simulated spectra are scaled to match. 

4. Results 

Using our scaling relations, we estimate the error on the values of T_i2 from our 
fiducial simulation using the uncertainties on our fiducial parameter values taken from 
the literature. The error budget is listed in Table [3 We find the largest contribution 
is from the uncertainty in the gas temperature. There is also a large contribution from 
the uncertainty in the effective optical depth, especially at z = 2. Interestingly, the 
uncertainty in f2 m also gives a substantial contribution to the total error budget. This is 
primarily due to the sensitivity of T_i2 on f2 m due to changes in the r.m.s. fluctuations 
of the gas density at the Jeans scale. Other uncertainties are less important, in particular 
flbh 2 , o% and h. Using our fiducial model, we find the Lya effective optical depth of the 
IGM at z = [2,3,4] is reproduced by T_i 2 = [1.29±g;|g, 0.86±g;||, 0.97±g;||]. 

We compare our results with other observational estimates of T_i2 in fig. El The filled 
triangles show the metagalactic hydrogen ionization rate computed from the updated 
UV background model of Madau, Haardt & Rees (1999), hereafter MHR99 (Francesco 
Haardt, private communication), based on the contribution from QSOs and re-processing 
by the IGM (e.g. HM96). Our data appear to be inconsistent with the IGM being kept 
ionized by QSOs at z ^ 4. For comparison, the filled squares give T_i2 calculated using 
the QSO rate above, plus an additional source in the form of Young Star Forming Galaxies 
(YSFGs). These ionization rates are in good agreement with our results, suggesting that 
a substantial contribution from galaxies appears to be required at all redshifts. 

Constraints on T_i2 from the proximity effect measurements of Scott et al. (2000) 
are shown by the hatched box. Our results are consistent with the lower end of these 
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Table 2. Percentage error budget for T_i2 from estimates of various cosmological and astro- 
physical parameters, listed approximately in order of importance. The total error is obtained by 
adding the individual errors in quadrature. 
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Figure 3. Comparison of our best estimate for T-i2 at z = [2,3,4] with constraints from 
observations. Our data are plotted with filled diamonds, and the grey shaded area shows the error 
bounds. The filled squares and triangles show the estimated contribution to the metagalactic 
ionization rate from QSOs±galaxies and QSOs alone, based on estimates from the updated 
model of Madau, Haardt & Rees (1999) including UV photons from re-processing by the IGM. 
The hatched box gives the constraint on T_i2 from the proximity effect (Scott et al. 2000) and 
the estimate from Lyman-break galaxies assuming a global spectral index of a = 1.8 is plotted 
with a filled circle (Steidel et al. 2001). The data point has been offset from z — 3 for clarity, 
and the upper error limit is not shown. The star shows the best estimate of r_i2 at z = 1.9 
from Tytler et al. (2004), including their errors from the uncertainty in r e ff, as and fib- 



estimates. The filled circle gives the T_i 2 we calculate assuming a spectral index of 
a = 1.8 using the metagalactic ionizing radiation intensity inferred from Lyman-break 
galaxies (Steidel et al. 2001). Finally we also plot the value of T_i 2 inferred by Tytler et 
al. in their study of the Lya forest at z = 1.9. The error bars include the uncertainties in 
r e ff, us and fib; based on scaling from their hydrodynamical simulations. Our data is in 
reasonable agreement with the Tytler et al. result, allowing for differences between the 
exact gas temperatures and numerical method. 
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5. Conclusions 

In recent years, compelling observational evidence has led to the acceptance of a stan- 
dard cosmological model which is flat, has low matter density and a substantial contribu- 
tion of vacuum energy to the total energy density. Within this 'concordance' cosmological 
model the current generation of hydrodynamical simulations predict values for the meta- 
galactic hydrogen ionization rate, required to reproduce the effective Lya optical depth 
of the IGM in the range z = 2 — 4, which are about a factor of four larger than those in 
an Einstein-de Sitter model with the same r.m.s. density fluctuation amplitude a%. The 
ionization rates estimated from the Lya forest opacity are more than a factor two larger 
than estimates from the integrated flux of optically/UV bright observed QSOs alone. 
This discrepancy increases with increasing redshift. We confirm the findings of Tytler et 
al. (2004) at z ~ 1.9 that the estimated ionization rates from simulations of a ACDM con- 
cordance model are in reasonable agreement with the estimates of the integrated ionizing 
flux from observed QSOs plus a significant contribution from galaxies as in the model of 
MHR99. The estimates of the ionization rate are also in agreement with the lower end of 
the range of values from the proximity effect. This new agreement appears to be in part 
due to the currently favoured low value of O m , which increases the required value of 
Fhi , and more accurate measurements of the Lya forest opacity. We also find that the 
uncertainty on the magnitude of Thi inferred from simulations is greater than previously 
estimated, with the primary contribution to this error coming from the temperature of 
the low density gas in the IGM. Better constraints on the thermal history of the IGM 
are required. Additional physics such as radiative transfer, galactic feedback and metal 
enrichment may need to be incorporated into simulations in a more realistic fashion. 
However, our estimate of the error on T_i2 most likely accounts for the modest changes 
expected from these processes, suggesting we have obtained a consistent constraint on 
the metagalactic hydrogen ionization rate. 
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